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Abstract: We study the strong coupling behaviour of fixed length single trace operators 
in the scalar SU(2) sector of = 4 SYM. We assume the recently proposed connection 
with a twisted half-filled Hubbard model. By explicit direct diagonalization of operators 
with length L = 4, 6, 8 we study the full perturbative multiplet of those lattice states which 
have a clear correspondence with gauge theory composite operators. For this multiplet, 
we follow the weak-strong coupling flow to free fermion states and identify in particular 
the precise asymptotic fermion configuration. Next, we analyze the Lieb-Wu equations of 
the twisted Hubbard model. For the antiferromagnetic state we derive its strong coupling 
expansion working at L up to 32. We also study the lightest state in the perturbative 
multiplet. This state is non trivial since involves complex solutions of the Lieb-Wu equa- 
tions. It is particularly interesting for AdS^ x duality since it is dual to the folded 
string semiclassical solution in the thermodynamical limit. We are able to perform the full 
analysis and compute the next-to-next-to leading terms in the strong coupling expansion 
for the non trivial lengths L = 12 and L = 20. A general formula is proposed for the NLO 
expansion for any L = A{2k + 1), A; G N. 
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1. Introduction 

The quantum behavior of = 4 super Yang-Mills in the planar limit is crucial in the 
context of AdS / CFT correspondence |l], ^, ^, |5| . In particular, the anomalous dimensions 
of certain single trace operators in the planar limit of the = 4 theory can be compared 
with the masses of string states on AdS^ x For instance, the comparison turns out 

to be particularly favorable for BMN states where the gauge-string matching can be 
done at the perturbative level in the planar gauge theory at the price of analyzing long 
composite operators with a large number of constituent fields and few impurities. 
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Apart from the gauge-string connections, the = 4 quantum theory has a rich internal 
structure suggesting its quantum integrabihty. The calculation of anomalous dimensions 
in specific sectors of the = 4 theory can be cast in algebraic form by computing the loop 
corrected dilatation operator Q . The huge mixing problem is then reduced to the analysis 
of the eigensystem of the finite dimensional matrix representing the dilatation operator. 

Remarkably, at one-loop, the dilatation operator can be identified with the Hamiltonian 
of the integrable XXX spin 1/2 lattice model |^]. In the SU(2) sector, Beisert, Dipple 
and Staudaucher (BDS) |l^ proposed a Bethe Ansatz for the 2-body 5-matrix in 

agreement with the explicit three loop dilatation operator and consistent with all loop BMN 
scaling [0]. The BDS equations describe a spin model of the Heisenberg-type with long 
range couplings. The range of the spin interaction increases as the loop order is increased. 
After having built the five loop BDS Hamiltonian, they could match the gauge theory 
predictions known up to three loops. A disagreement with the gauge theory calculation for 
operators with classical dimension L is expected precisely at L loop order due to wrapping 
terms. In the thermodynamical limit, these terms are negligible and an all-loop Bethe 
Ansatz was proposed. 

Using the BDS equations it is possible to compute the largest energy state on the 
chain In the Heisenberg model language, it the non trivial antiferromagnetic 

state. Its energy in the thermodynamical limit perfectly agrees with the ground state 
energy of a (twisted) one-dimensional Hubbard model at half-filling The connection 
between the spin model and the itinerant fermion model can be understood as follows. 
The spin model is nothing but the strong coupling expansion of the fermion model, where 
strong means the the hopping term is treated perturbatively. The effective Hamiltonian 
for the strong coupling of the Hubbard model contains interactions with a range increasing 
with the order of the expansion. This is because successive applications of the hopping 
operators connect lattice sites with increasing distance. 

With an impressive breakthrough, Rej, Serban and Staudacher (RSS) proposed the 
Hubbard model as the correct microscopic model behind the integrable structure of the 
M = 4 SYM dilatation operator [^]. In other words, they suggest that it could predict at 
all loops and non perturbatively the anomalous dimensions of the gauge theory operators 
for any finite L. This proposal also overcomes the problems related to the wrapping inter- 
actions in the long range spin model |15|. The RSS proposal is still a conjecture although 
with robust theoretical motivations. To pursue its assessment, it would be necessary to 
perform a four loop calculation in the gauge theory. Waiting for this check, the Hubbard 
model can be considered as a powerful Ansatz for the description of the gauge theory at 
finite operator lengths and beyond the perturbative regime. 

Actually, the appearance of the Hubbard model remains somewhat mysterious and 
intriguing |jl^ . In particular, the Hubbard model describes fermions with spin and admits 
states with two fermions in the same site. This extra states cannot be identified immediately 
with gauge theory operators. Indeed, it is not clear what is the role of its extra states 
which at strong coupling definitely mix with the perturbative states. Waiting for a better 
understanding of the role of Hubbard model, we can assume an optimistic attitude and 
exploit it to investigate the weak-strong coupling behavior of anomalous dimensions. 
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The main technical tool in the analysis of states in the Hubbard model are its Lieb- 
Wu equations which encode integrability. The Lieb-Wu equations are precious in the study 
of the thermodynamical limit including finite size corrections [p^ , [T^ , p!s[| . However, the 
machinery works well only for particular states. In the thermodynamical limit, the Bethe 
roots accumulate on a discrete set of non trivial curves in the complex plane. The integral 
equations for the root density are quite difficult due to unknown shape (and number) of 
the contours. At weak coupling, the Lieb-Wu equations reduce to the EDS Ansatz of 
Heisenberg type. In some favorable cases a solution for the thermodynamical limit can 
be found. Remarkably, the solution can also be matched with specific semiclassical string 
states |2^, 23, 24]. The analysis of the Lieb-Wu equations is more problematic and 

the full spectrum of the twisted Hubbard model at finite large L remains a quite difficult 
task. 

As an alternative approach to solving the Lieb-Wu equations, J. Minahan has recently 
proposed |25] to analyze directly the twisted Hubbard model Hamiltonian on small lattices 
to understand the features of the spectrum. This is an interesting approach aimed at un- 
derstanding the strong coupling behavior of states associated to gauge invariant operators. 

In this paper, we analyze single trace cyclic operators with zero SU(2) spin of the form 

Tr (Z'^ ^>-^) + . . . (1.1) 

with large J . Here, Z, and <I> are the charged scalar fields in the = 4 supermultiplet. As 
usual, dots stand for various other orderings of the scalar fields required to obtain eigen- 
states of the dilatation operator. The above set of operators, all with classical dimension 
2J is closed under perturbative renormalization. It is not clear what happens at strong 
coupling. Indeed, it has been suggested that some surprise could occur [^]. In this di- 
rection, it seems to be important to understand the mutual role of the perturbative states 
and the additional states in the Hubbard model. 

As a first step of our analysis, we collect and explain some general features of the states 
belonging to what we call the perturbative multiplet. These are the states in the Hubbard 
model with a clear identification with single trace operators in the gauge theory. Then, 
we exploit the direct approach of [|^] to study the full spectrum at L = 4, 6, 8. The case 
L = 4 is discussed in [^] and here it is reviewed to present the method and prepare for 
larger L. The cases L = 6 and 8 are more involved and reveal interesting features. In all 
cases, we determine precisely the fiow to large g. 

The direct analysis can hardly be pushed to much larger values of L. Hence, we revert 
to the numerical exploration of the complete Lieb-Wu equations. We perform the analysis 
of the antiferromagnetic state (the one with highest anomalous dimension) working with 
up to L = 32 sites and providing the NNLO strong coupling expansion of its anomalous 
dimension. 

Next, we consider the lightest state in the perturbative multiplet. Exploiting some 
features of this state in the thermodynamical limit, we solve numerically the Lieb-Wu 
equations at L = 12, 20 (the first non trivial cases of L = 8A; + 4, A; G N). These are 
remarkable values of the composite operators length. Nevertheless, the analysis is per- 
formed in full details obtaining also in this case the NNLO strong coupling expansion of 
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the anomalous dimension. This analysis is important in our opinion since it is an explicit 
study of the Lieb-Wu equations at large but finite L. 

The paper is closed by a conjecture about the general form of the lightest state for 
lengths L of the form 4(2A; + 1), A; G N. A discussion of the comparison with the BMN 
limit is also discussed. 



2. Anomalous dimensions in J\f — A SYM and the Hubbard model 

The four dimensional J\f = 4 SYM theory is finite. Its non trivial quantum properties 

are encoded in the behavior of gauge invariant composite operators. The gauge group is 
SU{N) and we arc interested in the planar limit N — >■ oo. We introduce the coupling g 
related to the large A'^ 't Hooft coupling A 

5' = ^, ^ = 9ymN. (2.1) 

In the so-called SU(2) sector of = 4 SYM, we consider gauge invariant composite oper- 
ators of definite scaling dimension 

TV {Z^-^ + . . . (2.2) 

where Z, $ are charged scalar fields. The classical dimension is L, the number of fields. 
In the following, we shall simply call L the length due to the lattice representation that we 
are going to introduce. The above composite operators have non trivial renormalization 
properties and acquire anomalous dimensions with all-loop corrections 

A{g) = L + Y.^i9''- (2.3) 

As we discussed in the Introduction, they can be computed as the eigenvalues of a charge 
D. It is the dilatation operator and belongs to an infinite tower of commuting charges. Its 
perturbative expansion (at two loops) is 

2 4 
V = L+^ ^(1 - ai ■ ai+i) - ^ ^(3 - 4ai ■ ai+i + ai ■ ai+2) + ■■■, (2.4) 

i i 

where the matrices act on the cyclic states of a spin 1/2 chain. A particular operator is 
mapped to a spin chain state in the natural way, e.g. 

Tr{ZZ^Z ■■■) — >\ TTiT • • •) + cyclic translations (2.5) 

The dilatation operator at finite g is thus non local. The RSS proposal identifies the up and 
down spin with fermions in two spin states. It also add states with two fermions occupying 
the same site. On this enlarged state space, RSS define an Hubbard-type Hamiltonian. In 
the following we shall consider the case L € 2N and M = L/2 and look at operators with 
zero SU(2) spin. The explicit local Hubbard Hamiltonian is 

H = Ho + -^Hi, (2.6) 
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where 

L 

Ho = L- ^n^,i ni^i, n^^ = c^^, (2.7) 

(2.8) 



i=l 



cr=T,i \j=l / 

The fermion creation and annihilation operators satisfy canonical anticommutation rela- 
tions 

{c.,.,c.'j} = .} = 0. (2.10) 

Periodic boundary conditions are understood. The twisting phase in the boundary link is 
fixed at </> = 7r/2. 

The Hamiltonian (p.6|) is symmetric under the SU(2) generators 



= ^EKi- M> (2-11) 

j=i 

^+ = E4.^i.' S- = {S^)^. (2.12) 



i=l 



At half-filling and with an equal number of up and down fermions the z-component is 
automatically zero. 



The Hamiltonian ( |2.6D is also invariant under the shift 

^ e^<^/^ c,j+i, j = 1, . . . , L - 1, (2.13) 
Ca,L ^ e^'^(^+i^/^c,,i, (2.14) 

with related transformation properties of the operators. The states invariant under this 
symmetry will be called cyclic states and are the relevant ones to represent single trace 
operators in the gauge theory. 

The Hamiltonian can be written in a simpler form after the transformation 

c, = e*"^Ac,. (2.15) 
The on-site Coulombian term Hq is unchanged. The hopping part Hi becomes 

L 

and the shift symmetry is written simply 

Cn^e^'^"^Cn+l. (2.17) 
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The hopping term is diagonahzed by introducing fermion operators in momentum space 

L 



where the lattice momenta can take the values 



Pn = ^, n = 0,l,...,L-l. (2.19) 
The dispersion relation for the scaled Hamiltonian Hi/\/^ is then 

r- ( 2-1171 IT \ , ^ 

e„ = y2cosf — + — j . (2.20) 

Cyclic states of the original Hamiltonian with L fermions can be built in momentum spaces 
by acting on the vacuum with L ot.p operators with a total momentum ^ p being an odd 
multiple of tt due to the phase factor e^^'^^^ . 

3. The perturbative multiplet 

As we have seen, the RSS construction introduces additional states with fermion double 
occupancy which do not have a direct correspondence with the gauge theory composite 
operators. As we remarked in the Introduction, the role of these states is not totally clear. 
However, we can exploit them as an auxiliary device and focus on what we shall denote 
the perturbative multiplet. This is the set of states which at g — > reduce to states with 
no double occupancy and have A(g) — > L, the maximal value at g = 0. In other words, 
these are the states which flow at weak coupling to states that can be naturally mapped 
to gauge invariant single trace operators. As g increases, they are no more the maximal 
energy states and mix with all the extra states. 

At very large g, any state flows to a free fermion state which is an eigenstate of 
the hopping term Hi. Our aim is to understand the asymptotic free fermion content of 
specific states in the perturbative multiplet. In particular, at large g we find A{g) ~ 
S g + 0(1) where 6 is the hopping energy of the asymptotic free fermion state. We show 
in Fig. (0) the qualitative description of the spectrum. The up most state is the so-called 
antiferromagnetic state (AF). The bottom part of the perturbative multiplet is composed 
of what we shall call light states, where light means that the anomalous dimension is small 
in the weak coupling region. The multiplets of extra states with double occupation are also 
schematically shown. In principle they can cross the perturbative multiplet. 

In the following sections, we shall discuss some general features of the AF and light 
states that can be derived from general principles and an analysis of the BDS equations. 

3.1 The antiferromagnetic state 

At even L the AF state is non degenerate. Its anomalous dimension is known in the 



thermodynamical limit. In terms of the planar coupling A = Svr^ g'^, it reads |13] 



lim 
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where 



This function is well known to be non-analytic in x = 0. However, it admits the asymptotic 



expansion |27| 



1 ^ 

m = i - ^ + ^ ^,^x'"^ + (3 3) 

TT 4 ^ — ' 

m=l 

where 

(2m - l)(22"^+i - 1) [(2m - 3)!!]=^ C(2m + 1) , ^,„_^ ,„ 
= 23>n-i(m-l)! = (3-4) 

Despite being only asymptotic and not convergent, the above expansion has been shown to 
reproduce correctly the second order perturbative correction at large A. This means that 
we can expand the anomalous dimension at large A as 

Aaf(A, L) = ai{L) VA + a2{L) + a3(L)-^ + . . . , (3.5) 

and the limits limi_^oo CLk{L) are obtained by replacing /(x) by its asymptotic expansion. 
At second order, we have 

/(x) = --f + + (3.6) 

vr 4 4 VT'^ 

and we obtain 

1- KTH^L) ,^^,fn\ VA 3 7C(3) 1 

= ^^ + Z + ^C(3)i + .... (3.7) 
vr 4 87r3v2 g 



Here, A^p°'^^(A,L) stands for the expansion Eq. (3.5). Later, we shall compare this pre- 



diction with the finite L analysis of the Lieb-Wu equations. 
3.2 The light states 

At the bottom of the perturbative multiplet there are light states. These are light in the 
sense that their anomalous dimension is small in the perturbative region. These states are 
highly non trivial in the BDS description. In the Heisenberg language, up to a change of 
sign in the anomalous dimension, they are states which can be built adding many excitation 
over the antiferromagnetic state respecting the constraint of zero spin and cyclicity. Not 
very much is known about these states since they correspond to non trivial distributions of 



Bethe roots in the thermodynamical limit |28|. In some cases, the Bethe Ansatz equations 



can be solved at L ^ cxd. An example is the lightest state which is associated with a 



limiting double contour distribution! 20, 21]. For brevity, we shall denote this state as |FS) 



since in the BMN limit it is dual to the so-called folded string solution |2£] 
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The state |FS) can be studied at finite L and its anomalous dimension can be loop 
expanded. If we express the result in terms of the 't Hooft coupling A we find 

Afs(A,L) = L + J;^A^. (3.8) 

The reason for the various explicit powers of L is that all the coefficients ci{L) have finite 
limits as L — > oo. Hence, the state |FS) admits the BMN limit 

L ^ oo, fixed = 87r2 (^^^ . (3.9) 

It is usual to introduce the coupling A' = A/J^ where J = L/2 is the angular momentum 
of the folded string. The coupling A' is fixed in the BMN limit. The BMN limit of the 
anomalous dimension is then 

lim \ Afs iX'L'^,L) = F(A'), A' fixed. (3.10) 



The function -F(A') has been first computed in |^. At small A', it reproduces the gauge 
theory perturbative expansion 

^,y) = i + ^A' + .... (3.n) 

At large A', F(A') ~ (A')^/^, the typical behavior expected from AdS/CFT duality |0|. 

In this paper, we work at finite L and cannot access the limit Eq. ( |3.1[1| ). Instead 
we are studying the |FS) state at fixed L, expanding Afs at large A'. This is the same 
procedure we followed for the AF state and is the kind of investigation described in |p5|| . 
We remark that there can be important differences between this and the BMN limits. At 
finite L and A' there can be terms with ambiguous A', L — > oo limit. An example could be, 
for instance, 

^'/L (3 12) 

l + A'/L' ^'^■^^^ 

tending to 1 when A' — > oo and to when L ^ oo. 

BMN scaling appears to be a general feature of the light states and is quite effective 
in the search for the corresponding dual string states. It is natural to look at BMN scaling 
as an infinite volume limit where the fixed ratio -v/V ^ g/L is interpreted as a finite size 
scaling variable. This smooth thermodynamical limit on the lattice can be explored more 
explicitly by evaluating in the Heisenberg model the SU(2) invariant correlation function 

Gk = {FS\ai-ai+k\FS). (3.13) 

In Fig. (^) we show its behavior at various L. Details of the calculation are reported in 
App. A. The correlation function expressed in terms of the scaled position k/L — 1/2 
tends to a smooth curve at large L. Similar results can be obtained for the other low lying 
states in the perturbative multiplet. For completeness, we also show in Fig. (^) the spectra 
of one loop anomalous dimensions for all 5 = cyclic states of the Heisenberg model with 
even 8 < L < 16. 
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4. Analysis of the full spectrum 



In this Section, we begin the analysis of the information that can be derived in the frame- 
work of the full Hubbard model. In order to understand the general features of the weak to 
strong coupling flow we analyze the full spectrum at variable g and L = A,Q, 8. We follow 
the direct approach of Minahan [^] already applied to the case L = 4 that we also review 
to fix the approach, extending it to somewhat larger values of L. Later we shall discuss a 
different approach based on the numerical solution of the Lieb-Wu equations. 

4.1 Review of the L = 4 case 



As explained in the very nice investigation |25], in the L = 4 case and after restricting to 



cyclic states with 5 = 0, there are only 6 remaining states. The antiferromagnetic state 
is non degenerate. Its perturbative expansion involves only even powers of g as an exact 
discrete symmetry of the model 

A = 4 + J]q5'^. (4.1) 

The full spectrum can be easily evaluated numerically and leads to the weak-strong coupling 
flow shown in Fig. (^. 

As explained in [^5| , the perturbative expansion of A can be recovered quite efficiently 
from the expansion of the secular determinant 

P(A, g) = det (^Ho + ^ - A^ . (4.2) 

One finds, 

P(A, c/) = A^ - 17 A^ + (119 - 16c/2) A^ + (-439 + ITGc/^) A^ + (900 - 7165^ + 32/) A^ + 

+ (-972 + 1276^2 - 160/)A + 432 - 8i0g^ + 200/. (4.3) 

Replacing the expansion ( |4.1| ) and matching the coefficients we find immediately 

A = 4 6 5^ - 12 / 42 / - 318 / 4524 g^^ - 63786 g^"^ + 783924 - 8728086 g^^ + 

+ 93893622 g^^ - 1038217494 + 12181236666 g^^ + ■ ■ ■ (4.4) 

At strong coupling (large g) the leading behavior of the eigenvalues is linear in g with 
a slope given by the eigenvalues of Hi/^/2. These are not immediately obtained from 
the dispersion relation because not all multifermion states are allowed by the cyclic and 
5 = conditions. The explicit eigenvalues of the 6x6 matrix i?i/\/2 can be computed 
analytically and are 



0, 0, -2 Y 2 - \/2, 2 y 2 - V2, -2 V 2 + V2, 2 y 2 \/2. (4.5) 
From the free fermion dispersion relation Eq. ( p. 20 ) we see that the highest state has a 



slope that unambiguously identifies it with the free fermion state with the following level 
occupation 

nt = (0,3), ni = (0,3). (4.6) 

The notation means that the components of rio- are the energy modes of the a type fermions 
according to Eq. ( 2.2C| ). This is nothing but the ground state of the Hubbard hopping term. 



-9- 



4.2 Extension to longer operators 

The extension of the above direct approach to longer operators is in principle straightfor- 
ward. However, some technical issues must be clarified in order to make the procedure 
systematic. We now illustrate the general features and then discuss the L = 6 and L = 8 
cases. 

To build the relevant states we first enforce the S = condition. Since N-^ = N^, we 
have automatically = 0. The operators are ineffective on doubly occupied sites. 
Also, they cannot move around the unpaired fcrmions. As a consequence, we can partition 
the problem of listing S* = states according to several sectors where we fix (a) the number 
and positions of the paired fermions, (b) the positions of the unpaired fermions. 

Given such a sector, the fermions are no more itinerant from the point of view of the 
spin calculation. Then, spin zero states are simply the spin zero component of the SU(2) 
decomposition of the product of fundamental representations. The actual wave function 
of these states is obtained by taking independent antisymmetrizations with respect to pairs 
of up and down fermions. The independent antisymmetrizations are explicitly given by the 
SU(2) Young tableaux with two rows and A^/2 columns. The Young tableaux entries in 
each column determine the independent antysymmetrizations. 

To give an example. Suppose that we have 3 + 3 fermions on a L = 6 lattice. In a 
sector where there is one paired couple in the rightmost site and the unpaired fermions are 
in the 1, 2, 3, 4 position, we have 6 states 



l*,*,*,*,o,Ti), 



* =T or i, 



(4.7) 



with 3+3 fermions in total. Then, the relevant two Young Tableaux with their associated 
antisymmetrization prescriptions are 



1 


2 


CO 


4 



1 


CO 


2 


4 



(4.8) 

giving two spin zero states, in this sector. 

The number of Young tableaux with 5 = constructed with L/2 up and L/2 down 
fermions is the first coefficient in (as usual (— n)! = when n G N) 

i8i2Ar 



= e 

s=0,l,2,.. 



(2s + 1)(2A)! 



{N + s + 1)\{N - s) 



(4.9) 



where [s] is the spin s representation of SU(2). Hence the desired number is the Catalan 



number C 



L/2 



Cl/2 « 



Cr 



(4.10) 



'■^/^ (L/2)!(L/2 + l)!' 
Summing over sectors with p doubly occupied sites, we find the total number of S = 




{L-2p)\ 



L!(L + 1)! 



L-p \ 

L-2pj (L/2 - p)!(L/2 -p + l)\ ((L/2)!(L/2 + 1)!)2 ' 



(4.11) 
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This number is reduced roughly by the factor 1/L after the cychc projection. This is im- 
plemented rather easily as follows. We denote by U the unitary operator which implements 
the shift symmetry. It can be checked that = 1 on half- filled states. Cyclic states \s) 
satisfy U\s) = \s). We then consider the shift symmetrizer 

5 = i(l + [/ + ••• + [7^-2 + f/^-i). (4.12) 
Ij 

If TLq is the space of zero spin states and Co the space of cyclic zero spin states, it is clear that 
the image STLq contains a basis of Cq. Indeed, any state \s) G Cq can be written |s) = S\s) 
and thus belongs to is in STLq. Hence, we can compute the image STLq and apply the 
Gram-Schmidt orthonormalization algorithm to simultaneously produce an orthonormal 
basis and also remove linearly dependent states. Of course, this procedure can be applied 
to the states in each of the sectors that have been identified in the construction of S" = 
states. This means, we repeat, sectors with fixed paired fermions and fixed positions of 
unpaired fermions. In addition, the cyclic structure of the fermion configurations in the 
cyclic states greatly helps in performing the orthonormalization by restricting to states 
with the same configurations modulo translations. 

After these technical remarks, we analyze in turn the L = 6, 8 cases which, as we shall 
discuss, illustrate in our opinion some interesting features valid in more complicated cases. 

4.3 The non degenerate L = 6 case 

After the spin zero and cyclic constraints, there is 1 state with no paired fermions. Thus, 
the antiferromagnetic state is again non degenerate. Also, there are 10 states with two 
pairs, 14 states with three pairs and 4 states with all fermions paired. The total dimension 
is 29. Fig. (^) illustrates the weak-strong coupling flow of the spectrum. 

As compared with the L = 4 cases, we observe several crossings of the coupling depen- 
dent levels. Such crossings are well known in integrable models, where they are explained 



in terms of additional conserved charges commuting with the Hamiltonian [31|. 

It is not feasible to evaluate the secular determinant P(A, g) at least if we do not want 
to resort to numerical evaluations. Instead, we can determine the analytical perturbative 
expansion of the highest eigenvalue by standard perturbation theory of non degenerate 
eigenvalues. Let -^q be the normalized eigenvector of Hq associated with the non degenerate 
eigenvalue L. Then, we set £q = L and iterate for n>l 

n—l 



V, 



n 



-Hi Vn-l + ^ £n-k Ipk, (4.13) 



k=l 



£n = -{lpO,Vn), (4.14) 

Ipn = 77 T (^" ^0 + Vn) ■ (4.15) 

Ho — L 

The last equation is evaluated in the subspace (V'OiV'n) = where the (pseudo) inverse 
operator (Hq — Eq)^^ exists. The perturbative expansion of the eigenvalue of Hq + "^-f^i 
is then 

s{g) = Y.en(^Y . (4.16) 

n>0 VV^y 
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This algorithm is quite fast since it is based on matrix vector multiplications only. Applying 
it, we find 

2 o 4^63(7^ 621/ 7047^10 100953 g^^ , 

A = 6 + 6g —9q -\ h 

^^2 4 8 16 

2006127 g^^ 46992069 1100850183 g^^ 24465145473 g^o 

H \ h (4.17) 

32 64 128 256 ^ ^ 

514257122079 g^^ 10323764001117 g'^'^ 
-I 2 2 \ 

512 1024 



This expansion can be compared with the BDS approach based on the Heisenberg model |11]. 
The agreement is perfect up to five loop order which is where the long-range Heisenberg 
Hamiltonian is reliable at L = 6. As we remarked, the Hubbard model calculation is con- 
jectured to be exact at all orders in the loop expansions, although a proof is lacking. Of 
course, knowing the one-loop Bethe roots, the above expansion can also be obtained by 
perturbative expansion of the Lieb-Wu equations. We do not insist on this point, since we 
are mainly concerned with strong coupling properties. 

Again, we can exploit the direct diagonalization approach to identify the free fermion 
state to which the g = highest state flows. Comparing the slope of the highest eigenvalue 
(the maximum eigenvalue of i?i/\/2) with the dispersion relation Eq. ( 2.20 ) we find the 
two possibilities 

= (0, n, 5), = (0, n, 5). (4.18) 

where {n,n') can be (1,4) or (4,1). The contribution of these two fermions cancels in 
the energy. Indeed, ei + £4 = 0. This means that we flow to an excited state of the full 
Hubbard model. This is a clear consequence of the cyclic projection. Indeed, the ground 
state of the Hubbard hopping term for L = 6 is not cyclic and is instead odd under the 
transformation Eq. ( |2.13D . 

The spin zero condition determines uniquely the correct combination of states which 
is the antisymmetric combination 

-^(|0,l,5)t |0,4,5)i - |0,4,5)t |0,l,5)i) (4.19) 

As a check, we see that the largest eigenvalue of the explicit 29 x 29 matrix Hi on the 
cyclic spin zero states is non degenerate. 

4.4 The degenerate L = 8 case 

On a -L = 8 lattice there are 4900 half-filled states in the full Hubbard model. After 
the spin zero and cyclic constraints, there are 3 states with no paired fermions, 35 states 
with two pairs, 108 states with three pairs, 70 states with four pairs, and 10 fully paired 
states. The total dimension is thus reduced to 226. This is a remarkable reduction, but the 
dimension remains rather high. Nevertheless, we shall be able to complete the analysis. 
The maximum eigenvalue of Hq is threefold degenerate. It contains the antiferromagnetic 
state and other two states with lower anomalous dimensions. We postpone the discussion 
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of flow. Again, it is clearly not feasible to evaluate the secular determinant. Also, we must 
deal with the complication that there are 3 states with eigenvalue L = 8 at g = 0. 

We can determine the perturbative expansion of the highest eigenvalue by quantum 
mechanical formulae for perturbation theory of degenerate eigenvalues. In the case under 
consideration the degeneration is removed at second order in g. A very simple practical 
algorithm is then the following. 

Let Pq be the projector onto the degenerate eigenspace Sq with eigenvalue Eq = L = 8. 

Let 

P = Pogi(l-Po) ^ {l-Po)HiPo (4.20) 

be a 3 X 3 matrix restricted to 8q. Let its three eigenvectors be V'O) V^o; V'o some order. 
They have distinct eigenvalues. We iterate for n > 

1 / ^""^ A 

V'2n+l = TT _ ^ £2n+l-k " Hi V'2n , (4.21) 



Hq — Eq 

2n 



,fc=l, odd 



V2n+2 = ^ £2n+2-k'>Pk - Hitp2n+l, (4.22) 
k=2, even 

e2n+2 = -(■0O,^^2n+2), (4.23) 

V'2n+2 = 77 FT (e2n+2 ^0 + V2n+2) + "2^+2 V'o + a2n+2 V'o • (4-24) 

Hq — 1^0 

The coefficients are fixed by the condition Pq f 2n-2 = 0. 

The explicit matrix P is rather complicated. Its eigenvalues are the three roots of the 
equation (of course in agreement with the one loop calculation in [^) 

+ 40A2 + 464A + 1600 = 0. (4.25) 

Finding numerically the eigenvectors and applying the above algorithm with the three 
possible choices for ^q, we immediately find the three perturbative expansions 

A = 8 + 11.3022 g^ - 22.1706 / + 79.5035 / - 352.94 / + (4.26) 
+ 1777.24 5^° - 9743.47 g^'^ + 56739.6 g^^ - 825617 + • • • , 

A' = 8 + 5.45222 g'^ - 7.94042 / + 31.2193 g^ - 159.093 / + (4.27) 
+ 896.064 g^° - 5378.33 g^'^ + 33796.9 g^^ - 222137 g^'^ + ■ ■ ■ , 

A" = 8 + 3.24559 5^ - 1.88899 / + 1.2772 / + 1.0331 / + (4.28) 
- 8.2996 + 27.3006 g^^ - 70.9279 g^^ + 159.235 g^^ + ■ ■ ■ 

They are unavoidably numeric since they involve the algebraic irrational A. 

Again, we can compare with the BDS prediction at five loop. It is given in [0] terms 
of an algebraic number (denoted ip) which agrees with the above three A roots. The 
agreement is complete. Of course, beyond five-loop, the above expansions derived in the 
Hubbard model are new and must be checked against gauge theory perturbation theory. 
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Coming to the analysis of the weak-strong couphng spectrum flow, we see that it 
is quite complicated as illustrated in Fig. (^). The identification of the asymptotic free 
fermion state is less easy but straightforward. Let us denote by 

|A), |A'), |A"), (4.29) 

the three states with eigenvalues expressed by the above expansion. The highest state |A) 
is the antiferromagnetic state and does not cross the other states along the flow. Instead, 
the other two degenerate states undergo several crossing as g is increased. However, it is 
easy to follow them along the crosses. We enumerate states starting from the highest. The 
two subleading states turns out to be 3rd and 8th asymptotic free fermion states. This is 
more clearly illustrated in Fig. (^) where we show the first 9 ordered eigenvalues. Several 
crossings can be observed and in the end, the asymptotic eigenvalues remain separated. 
Remarkably, at large g the states |A') and |A") are close (in energy) to partner states |A') 
and I A") that we shall now discuss. 

As a first step toward the identification of the asymptotic states with free fermion 
states, we analyze the eigenvalues of the matrix Hi/\/2 and compare them with the dis- 
persion relation. If we denote by s, s', s" the asymptotic slopes of the energies of the three 
states I A), |A'), |A"), with respect to the coupling g, we find the following (unique) match 
in terms of the free fermion energies Eq. ( |2.20D , 

s = 2(eo + ei +£6 + £?), (4.30) 
S' = 2(80+87), (4.31) 
s" = 2(81+87). (4.32) 

The first relation permits to conclude that the highest states is nothing but the Hubbard 
model hopping term ground state. The other two relations identify two levels which are 
occupied by an up-down doublet t| in momentum space. The remaining 4 fermions (two 
up and two down) must be placed in the remaining levels with a total zero additional 
energy and respecting the cyclicity and zero spin conditions. This can only be achieved by 
leaving the four fermions unpaired and placing them symmetrical around the zero energy 
value. The level population is shown in Fig. (|8|). 

The first state on the left is the Hubbard model hopping ground state, as discussed. In 
the other states we have shown the two levels which are unambiguously filled with a pair ||. 
We have also shown a particular admissible distribution of the remaining two up, and two 
down fermions on the allowed four symmetrical levels. We have circled them with dashed 
ellipses to emphasize that this is just one component of the exact state. Indeed there are 
several possible distributions of the unpaired fermions. To further analyze, we take into 
account the spin zero condition. For both |A') and |A") the allowed states reduced to the 
two independent states that are obtained by antisymmetrizing two pairs of up and down 
fermions. 

Of course, there are two states because of the SU(2) decomposition 

1111 , , 

-®-(8)-®- = 2©l©l©l©0©0. (4.33) 

At At ^ ^ 
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To be explicit, in the case of |A'), two orthonormal states can be taken to be 




2(1 1,I4,T,i,T,o,o)p-| |,|,i,T,T,i,o,o)p + 



(4.34) 



-I l4>T4,i,T,o,o)p + | |,|,T,i,T,i,o,o)p), 
e'2) = ^(2| I44,i,T,T,o,o)p-| |,|,i,T,i,T,o,o)p-| |,|,i,T,T,i,o,o)p + 

-I :,tT,i,i,T,o,o)p-i :,:,T,i,T,i,o,o)p+2| :,:,T,T,i,i,o,o)p) (4.35) 



where the states \ - ■ ■)p are labeled with the fermion occupancy in momentum space and 
the momentum sites are ordered from the largest En (^o) to the smallest {e^ = —Eq). With 
this notation, the Hubbard model hopping term ground state is | J, J, J, J, 0, 0, 0, 0). 

Introducing analogous orthonormal states |e'/2) for the sector spanned by |A") and 
I A"), the free fermion asymptotic states associated with |A') and |A") are suitable linear 
combinations of \e[ 2) and \e'(2) that can be determined by perturbation theory in Hq. In 
both cases, the relevant asymptotic state is the one with (slightly) larger energy as can be 
seen from Fig. (^. 

At first order, the degeneration is not removed. We find the same constant shift in 
both doublets. At second order, we find a non trivial energy splitting. We do not report 
the expression of the asymptotic eigenvectors which is really not useful. Instead, we give a 
closed form for the split eigenvalues. 

If we denote by A'_|_ and A'j. the strong coupling expansion at second order of the 
doublets eigenvalues, we find 



A'± = 25 (eo + £7) + y + - 5± + ■ • • 



(4.36) 



where 




(4.37) 



(4.38) 



The expansion of the second eigenvalue is instead 



Al = 2g{£i + e'r) + — + -5'l + 



(4.39) 



where 




(4.40) 




(4.41) 
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Notice that this is correct for g > 0. Indeed, in general, the above strong coupHng ex- 
pansions must be written with g — > \g\ to respect the exact g — > —g symmetry of the 
spectrum. 

As a final comment, we remark that the above complicated expressions have been 
checked explicitly against the numerical evaluation of the levels with full agreement. 

4.5 Summary of the results for L = 8 

In conclusion, taking the upper states and evaluating the energy levels, we find the following 
result in the case L = 8. The AF state has been already discussed. The other two states in 
the perturbative multiplet have the following strong coupling expression of the anomalous 
dimensions 

A'+ = ^2 (^4 + 2V2 + ^20 + 14v^ ^ + ^ + i . . . , (4.42) 
A'; = Afs = 2\j 2 + ^/2 + ^ ^ + ^ + 1 + . . . (4.43) 
5. Direct analysis of the Lieb-Wu equations 

A posteriori, we can make some general comments on the previous analysis based on direct 
diagonalization. We are considering an Hamiltonian of the form 

H = HA + gHB, (5.1) 

where 5 is a coupling. As g flows from to 00, each eigenstate of H flows from an eigenstate 
of Ha to an eigenstate of Hb- It is clear that exact diagonalization permits to follow the 
flow for generic H^, Hb, however the method is unrealistic for large dimension of the 
Hilbert space. The problem is very general and, as such, has no simple solution. Of course, 
what saves the day in our context is integrability. In principle, the Lieb-Wu equations can 
be solved for a particular state, i.e. typical configuration of Bethe roots, without the need 
for huge calculations of eigensystems. Extending the calculation from g = 0, where the 
Bethe roots are those of the Heisenberg model, up to large g should permit in principle to 
determine the strong coupling behavior of any state. 

However, it is also clear that the general task of solving the Lieb-Wu equations for all 
states in the perturbative multiplet and fixed L (possibly large) is not easy [^2|. However, 
there are exceptions. These are the states where some general knowledge is available about 
the limiting distribution of Bethe roots at large L. In the next Sections, we shall discuss 
two important examples. The first is the AF state. At half-filling, it is the unique state with 
completely real solutions to the Lieb-Wu equations. The second example is the lightest 
state |FS). Here, we know that at large L the Bethe roots condense on two symmetric 
curves in the complex plane and we can exploit this information to evaluate them at least 
numerically. This case is considerably more difficult than the AF state because the Bethe 
Ansatz solution is complex. 



-16- 



Notice that in principle, one could use the original all-loop BDS equations. This 
calculation would be reliable in the thermodynamical limit including finite size corrections 
if needed. In this paper, we are concerned with finite L properties, and therefore we have 
explored the explicit (numerical) solution of the more difficult Lieb-Wu equations. 

5.1 Real solution: the antiferromagnetic state 

At half filling, the AF state is described by the only genuine real solution of the Lieb-Wu 
equations. They read 

L/2 

L = 2tt In + tan"^ 2{uj -V2g sin(g„ + 4>)) , (5.2) 

L/2 L 

27r Jfe = 2 ^2 tan~^('Ufe — Uj) — 2 ^ tan~^ 2{uk — V2g sm{qm + (f))) , (5.3) 

j=l m=l 

where, in our problem, n = 1, . . . , L, and k = 1, . . . , L/2. We focus on the case N = 4p 
where (f) = 7r/(2L) and the Bethe quantum numbers are 

{/„} = {0,1,2,...,L-1}, (5.4) 

r 2p-l 2p-3 2p-3 2p-l\ 



The iterative solution of the above equations is quite stable, as it is common when dealing 
with real solutions. Following the evolution at 5 oc of the energy of the highest state, 
we have checked that it flows at strong coupling to the ground state \iIjq) of the Hubbard 
model hopping term. In momentum space, this is the state where all positive energy levels 
are doubly occupied 

iV'o) = n' n <pjo)- (5-6) 

n=l (T=T,i 

The AF state remains non degenerate at all couplings and we can apply first order pertur- 
bation theory to determine the first subleading correction at large g. Also, from the nu- 
merical solution of the Lieb-Wu equations, we can evaluate the finite L next-to-subleading 

correction ~ 1/g. 

Summing up, our result for the expansion of the anomalous dimension of the AF state 
at finite L and large g reads 

AAF(ff,i^) V2 ^3 1 

The first term is the energy of the Hubbard model hopping term ground state with twist. 
Its explicit formula comes from the sum 



n=-L/4 ^ / 2L 



V2 
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The coupling independent term is a universal constant. It is evaluated by computing the 
matrix element 

j{'^o\^a\pa^^pa\gai^q\'i{;o) = L- -l-j = ^L. (5.9) 

p,q ^ ^ 

The next correction takes the numerical values reported in Tab. (|l|). The above result 



L 8 


12 


16 


20 


24 


28 


32 


5l 0.0250979 


0.0245631 


0.0243433 


0.0242308 


0.0241652 


0.0241234 


0.024095 



Table 1: Coefficient of the second order correction to the energy of the AF state at finite L E AN 
and large g. A simple polynomial extrapolation at L — > oo gives limL^oo = 0.0240(1). 

is an exact expansion in inverse powers of g at fixed L. As a non trivial check, it can be 
compared with the result Eq. (|3.7|) by taking the L — > cxd limit term by term. We obtain 



L-+00 L vr 4 g 



This is in full agreement with Eq. ( |3.7| ) since 

7 

6oo = lim 5l = — C{3) = 0.0239866. (5.11) 
5.2 A complex solution: the |FS) state 

The general geometry and symmetry of the Bethe roots for the |FS) state at large L are of 
great utility in determining them at finite L £ 4N. As we shall discuss, there are numerical 
difficulties when L = Ap with even p, i.e. Lmod8 = 0. Instead the case L = Ap with odd 
p, i.e. Lmod8 = 4, can be treated successfully. In the next Sections we shall present our 
detailed results for the non trivial cases L = 12, 20. 

5.2.1 L = 12 

At L = 12 there are 14 states in the perturbative multiplet and many more extra states in 
the full Hubbard model. We have determined the one loop Bethe roots for the |FS) state 
as we now discuss. The six Bethe roots are non zero and symmetric under u — > —u. This 
reduces the problem to solving three polynomial equations in three variables. Four roots 
are expected to be complex and the other two real. Applying the resultant technique we 
find that the complex roots are among the solutions of the following polynomial 

R{x) = 1 - 291816x2 + 2476695728 - 4740875459840 + 2716015001869568 x^ 

-587934012140484608 x^° +53336517749102178304x^2 - 1893188143985026269184 x^^^ 
+27036104777708125093888 x^*^ - 119124909860572824600576 x^^ 
-345397582972412910108672x2° + 1499936486421645590790144x22 
+8744914427777415217414144x2^ - 1155658862325646445510656 x2'5 (5.12) 
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+7466630646963993812402176 x^*^ - 273168717774897644721668096 a;™ 
+46033164223912000241008640x^2 - 219403347416043302531629056x3^ 
+965423684859142279840923648 x^^ + 4267809644452908595622707200 x^^ 
+886394409006105612647399424 - 2540126101042720433221140480x^^2 

-10295374756958736948626718720 x''*^ - 27408343376808874282372300800 x*^^ 
-16677333147947189938421760000 x"^^ + 23395667776149943429890048000 x^° 
+48702367142746624075235328000 x^^ + 43027672720198395903344640000 x^^ 
+21488277238759969133690880000 x^*^ + 5377757322242364014592000000 x^^ 
+496408368206987447500800000 x*^" . 

This polynomial contains several spurious solutions and is perhaps not the most econom- 
ical choice. Nevertheless, it contains the exact roots and we quote it for the reader's 
convenience. ^ The roots associated with the |FS) state are 

ui = —U2 = a, (5.13) 

Us = — U4 = a, (5.14) 

W5 = -^^6 = /3, (5.15) 

where 

a = 0.6762450414055523 + 0.9936333912043784 i, (5.16) 

(3 = 0.6780174422473694, (5.17) 

in agreement with the results in [^]. We start from this solution plus the condition for 
the momenta q 

gn = ^(n-l), l<n<L, (5.18) 

which is valid at 5 = 0. Then, we increase g and determine step by step the new solution of 
the full Lieb-Wu equations. This procedure is numerically stable and permits to determine 
the energy flow as well as the change in the Bethe momenta and u variables. 

We show in Fig. (|^, the evolution of the Bethe Ansatz solution {ui} as g is increased 
up to 5( = 15. More interestingly, we show in Fig. (|lO|), the evolution of the momenta. The 
flow permits to derive the asymptotic occupancy in the free fermion limit. In terms of the 
indices n, the final occupation of states is as follows. There are singly occupied modes with 
mode indices 

n = 0,11,2,9,3,8,5,6, (5.19) 



and two doubly occupied levels at modes 



n = l,10. (5.20) 



single polynomial for the case L = 8 is quoted in However, there is a misprint in one coefficient. 
The correct resolvent is R{x) = -1 + 648 a;^ - 36464 + 81664 x^' - 16128 x^ + 460800 x^° + 552960 x^^ . 
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As in the L = 8 case, the singly occupied levels contain 4 up fermions and 4 down fermions 
in a 5 = combination. There are several possibilities and only one can be selected by 
perturbation theory. We do not pursue the strong coupling correction analytically. Instead 
we determine the leading term at large g and the subleading contribution that can be 
obtained by first order perturbation theory. The leading contribution (the coefficient of g 
at large g) is 



2(ei + eio) = 2V2 ^cos ^ + cos ^ j = 4 cos ^ = 2^ 2 + ^2 + \/3. (5.21) 

The subleading term can be computed analytically and is 26/3. Hence we have found that 
for the |FS) state at L = 12 we have 



Afs = 2^/2 + \/2 + V3 g + ^ + 0.597(1) - + •••, (5.22) 

6 g 

where we have also indicated the fitted coefficient of the NNLO term in the strong coupling 



expansion. The agreement with the calculated energy is shown in Fig. (11). 

The asymptotic fermion configuration is shown in Fig. (|l^) where we simply draw the 
single particle level and their occupation without specifying the spin of the singly occupied 
levels. 



5.2.2 L = 16 

We can repeat the analysis for L = 16. In this case, we failed to obtain an exact resultant 
encapsulating the exact one loop Bethe roots. Instead, we have computed them numerically. 
The symmetry of the 8 roots is 



Ui = —U2 = Oi, 


(5 


U3 = —U4 = a. 


(5 


U5 = -uq = 13, 


(5 


u-j = -Us = 13, 


(5 



where 

a = 0.9011983985707239 + 0.5000879064837407 i, (5.27) 
[3 = 0.915478863907937 + 1.4850185722704357 i. (5.28) 

The imaginary part of a is quite near to ^. This is a source of instability in the solution 
of the Lieb-Wu equations. Indeed, as g is increased, we numerically observe that four of 
the Bethe roots tend quickly to a singular configuration. This problem does not occur if 
L mod 8 = 4. We do not try to deal with the singularities of the L = 16 case and instead 
study the more involved, but more stable case L = 20. 
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5.2.3 L = 20 

The 10 one-loop Bethe roots satisfy the following conditions: 



ui = -U2 = ae M>o, (5.29) 

ua = —U4, = P, (5.30) 

U5 = -uq = (3, (5.31) 

m = -Us = 7, (5.32) 

ug = -uio = 7- (5.33) 

We find the following numerical solution 

a = 1.1309564538305164, (5.34) 
P = 1.1310261843923932 + 0.9998455911389437 i, (5.35) 
7 = 1.1784184821892867 + 1.980402937535511 i. (5.36) 

We show in Fig. (|l3|), the evolution of the Bethe Ansatz solution {ui} as g is increased up to 
g ~ 18. Fig. (^) shows the evolution of the momenta. Again, we can derive the asymptotic 
occupancy in the free fermion limit. In terms of the indices n, the singly occupied modes 
are 

n = 0, 1, 3, 4, 5, 6, 7, 9, 10, 11, 13, 14, 15, 16, 18, 19, (5.37) 
and there are again 2 doubly occupied levels at modes 

n = 2, 17. (5.38) 

The singly occupied levels contain 8 up fermions and 8 down fermions in a S* = combi- 
nation. The leading contribution to the anomalous dimension (the coefficient of g at large 
g) is 

^ / 9tt 697r\ vr 

2(62 + £17) = 2V2 (^cos — + COS — J = 4cos -. (5.39) 

The subleading term can be computed analytically and is 73/5. Hence, in summary, the 
|FS) state at L = 20 admits the strong coupling expansion 

Afs = 4 cos ^ 5 + ^ + 0.953(1) - + •••, (5.40) 
40 5 g 

where we have also indicated the fitted coefficient of the NNLO term in the strong coupling 
expansion. As before, we shown the agreement with the calculated energy is shown in 
Fig. (^. 

The asymptotic fermion configuration is completely similar to the L = 12 case. The 
doubly occupied levels are in the middle of the single particle positive energy levels. The 
other positive energy levels are singly occupied, as well as their mirror levels with negative 
energy. 
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5.2.4 Conjecture for the |FS) state at general L = 4(2A; + 1) 

The results at L = 12 and 20 are quite symmetric and completely similar. It is natural to 
conjecture that for all L = A{2k + 1), the pattern is identical. This means that the |FS) 
state is obtained at strong coupling as the state with the following properties. 

1. The positive energy single fermion levels are all occupied with one fermion, with the 
exception of the central levels with mode numbers 

n = k, L-k-1. (5.41) 

These are doubly occupied. 

2. The negative energy levels which are mirror of singly occupied levels are also singly 
occupied. 

3. The negative energy levels which are mirror of doubly occupied levels are empty. 

Evaluating the leading and subleading contributions to the anomalous dimension gives the 
strong coupling expansion 

Aps(5, L) = 4 cos ^ 5 + '"^^ + ' + O Q) . (5.42) 

Eq. ( ^.421) is expected to be the exact strong coupling expansion at any fixed L = 4(2A; + 1). 
We conclude with a comment about the other cases L = 4(2/c). The explicit results at L = 8 
and preliminary data at L = 16 suggest that the |FS) state is again flowing to a state like 
the above with the same expressions for the leading and subleading terms in the strong 
coupling expansion. Indeed, at L = 8, the above parametrization reproduces the exact 
result that we derived by exact diagonalization. However, we have not enough empirical 
support to firmly establish this result. 



6. Summary of results and discussion 

To summarize, we report our main results for the large g expansion at fixed L of the 
anomalous dimension of the antiferromagnetic and |FS) states. For the antiferromagnetic 
operator we have found 

J = . . n g + 7 + <^AF,L - + •••• (6.1) 

L L sm 2L 9 

The first two terms are exact. About the last one, we have shown how to compute it for 
large L. We have also provided the asymptotic limit 

For the folded string dual, we have found (at finite L mod 8 = 4) 

. r .^ . TT 3L2-2L + 8 , L , , 

Afs(5,^) =4cos — 5+ ^ + <5a,l - + •••, (6.3) 
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where we have exphcitly computed 



12 



0.0498(1), 6x,20 = 0.0477(1). 



(6.4) 



In particular, the leading and subleading terms of these expressions are exact at all finite 
L. As such they are beyond the region of applicability of the BDS approximation which, 
due to wrapping terms, is limited to L — > oo. Indeed, they are a genuine result provided 
by the Hubbard model framework where they express properties of specific free fermion 
states. 

As we remarked, our expansions are obtained by taking A' large at fixed L. An im- 
portant issue if then the comparison with the BMN limit. Indeed, as pointed out by 
Minahan |25|, it is not totally clear how to recover the (A')^/^ behavior of string states 
from the strong coupling expansion of the Hubbard model. The tricky proposal in |^5| is 
based on the assumption that at large g there are doubly occupied levels with small single 
particle energy of order 1/L. Unfortunately, we have seen that this is not valid for the 
considered states with many excitations. The doubly occupied levels give an asymptotic 
slope A/g which is of order 1 for the minimal energy state |FS). 

The BMN limit is obtained by fixing A' = A/J^ and taking J ^ oo where J = L/2 is 
the angular momentum of the semiclassical dual state. We have seen that at finite L there 
can be correction terms (like Eq. ( p.l2| )) with an ambiguity in the L, X' ^ oo limit. If we 
want to compare with the BMN limit we have to enforce the correct ordering and require 
A' <C L. Let us see the role of this constraint in the case of our data at L = 12 and 20. 
The BMN anomalous dimension of the folded string can be written 



lim 



AFs(Aa^L) 
L 



F(A') = ^Kiq) 



AqX' 



+ 



1 



1/2 



where q = q{X') is the solution of 



4^ 



7r2 {K{q) - E{q)f E{qf 



(6.5) 



(6.6) 



and K{q), E{q) are standard complete elliptic integrals of the first and second kind. 

We show in Fig. ( p!6[ ) the comparison of F(X') and the ratio A/L for the lightest state 
at L = 12 and 20. The left panel shows that there is good agreement for A' up to about 
3 — 4, reasonably within the BMN scaling window. This is a rather large value suggesting 
that the agreement is working beyond perturbation theory. Indeed, we show in the same 
figure the 8th and 9th order perturbative expansions of F{X') which read 



(6.7) 



v/\5 



F(A') = 

1 + 0.089004 A' - 0.013272 (A')^ + 0.002839 (A')^ - 0.000676 (A')'' + 0.000172 (A 
-0.000047 (A')^ + 0.000014 (A')^ - 4.315511 • 10~^ (A')^ + 1.424401 • 10"*^ (A')^ + . . . . 

The two curves suggest a convergence radius around A' ~ 3, somewhat smaller than the 
region of agreement. If so, this could be a signal that we are slowly recovering the BMN 
result. 
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7. Conclusions 



In this paper we have considered a particular class of gauge invariant operators in the planar 
limit of = 4 SYM. These are length L single trace operators in the SU(2) sector with 
zero spin. Remarkable members of this class are the so-called antiferromagnetic operator 
and the dual of the semiclassical folded string state. 

We have assumed a recently proposed relation between the gauge theory and a Hubbard- 
like model of itinerant fermions. This approach permits to evaluate the anomalous dimen- 
sions of the gauge invariant operators at all couplings. In particular, we access the strongly 
coupled region at fixed L. 

Our investigation has been based on two complementary techniques. First, we have 
evaluated the full spectrum of the model for operators with L = 4, 6, 8. This has provided 
useful information about the mutual relation between the states of the Hubbard model and 
the perturbative multiplet of states with a clean relation to gauge invariant operators. 

Then, we have investigated the numerical solutions of the Lieb-Wu equations. They 
are a powerful tool that permits, in principle, to follow a particular state from weak to 
strong coupling in a totally controlled way. Our results are very simple explicit formulas 
for the strong coupling expansion of the anomalous dimensions. They are expressed in 
terms of specific free fermion states whose properties are easily computable. 

We hope that this investigation will be useful to shed some additional light over the 
non perturbative features of the states in the multiplet as well as over the connection 
between the gauge theory and the underlying integrable Hubbard model. This work can 
be extended in several directions. Some are rather obvious like (a) 5 7^ states in the 
SU(2) sector, (b) sectors different than SU(2), (c) other particular states like for instance 
the dual of the semiclassical circular string solution. In principle, it should be possible to 
study the strong coupling region by direct perturbative expansion of the Lieb-Wu equations 
although this is a delicate analysis |Q. This could be a valuable effort. It would be very 
nice to reproduce in some limit of the Hubbard model the true strong coupling behavior of 
string states, i.e. the typical relation A ~ (A')"^/^ for the light states. This is non trivial at 
the numerical level since large A' and irrelevance of possible corrections like (|3.12| ) would 
require quite large lattice sizes L. 

The most interesting extension seems to be a detailed study of the other light states in 
the perturbative multiplet, perhaps exploiting in a deeper way their nearly-BPS natuie p^ . 
As L increases, our investigation shows that there is a growing number of states with a 
smooth L — > 00 limit. This is not unexpected and they should be described by a suitable 
effective theory in the continuum. This line of analysis have been discussed in [35| in the 
context of the loop-corrected Heisenberg model and should be extended to the Hubbard 
model |^6[. For these states, the finite L analysis of their associated Lieb-Wu solution could 
provide some hindsight on the possible limiting distribution of Bethe roots and suggest a 
strategy to evaluate their thermodynamical limit. In the end, this could lead to new 
examples of AdS/CFT specific dualities. Indeed, the search of new string states dual to 
novel Bethe Ansatz solutions seems to be far from the end |36]. 
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A. Exact diagonalization of the S = sector of the Heisenberg model 

In this Appendix, we report some useful techniques for exact diagonalization of Heisenberg- 
like models in sectors with fixed SU(2) spin. As discussed in the analysis of RSS, the two 
loop dilatation operator is 

V = L + — ^(1 - (Ti • (Ti+i) - — ^(3 - 4(Ti • fJi+1 + CFi ■ CFi+2) H , (A.l) 

i i 

where we assume periodic boundary conditions. This operator acts in the perturbative 
multiplets, i.e. on the 5 = cyclic states of the Heisenberg spin model. We can write the 
various ct, • Uj terms by means of transpositions operators flipping spins at sites i, j 

1 

The dilatation operator can be rewritten 



Pi^j = -{l + ai-aj). (A.2) 



P = L + ^2 ^ ( 1 - Pi,i+i ) + |- ^ (4Pi,i+i - - 3) + • • • (A.3) 

i i 

where we notice that (under periodic identification of the boundaries) 

Pi, i+2 = P'i+l, i+2 Pi, i+l , (A. 4) 

and the full operator is written in terms of elementary transpositions only, i.e. transposi- 
tions of nearest neighboring spins. Also, cyclic states are states invariant under a lattice 
shift T that can also be written in terms of elementary transpositions as the product 

T = Pl-i,l-- - P2,zPi,2. (A.5) 

We are interested in diagonalizing the dilatation operator in the S" = sector. The states 
in this sector are associated with standard SU(2) Young tableaux with two rows and L/2 
columns, as already discussed. However, now we are no more interested in the detailed 
spin positions and we do not need translating the YT in explicit (anti)symmetrized states. 
Instead, we can exploit an old computationally efficient parametrization of orthogonal 



states 1 37] which turns out to be very suitable for our problem. 

We associate a state |y) to each distinct spin zero Young tableaux Y. Then, the 
nearest-neighbor transposition P^^k+i has the following matrix elements 



(Y'\Pk.k+i\Y) 



a/ 1 — p^, y is obtained from y by fc -|- 1 (A-6) 
0, otherwise. 
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The number dk{Y) is the axial distance between the labels k and k + 1 inside Y. It is 
defined as the sum of steps which are required to move from k to k + 1. The steps are 
positive on the right and upward and negative otherwise. 

The above representation of states is not very convenient for the analysis of generic cor- 
relation functions. However, it is quite efficient for energies and SU(2) invariant correlation 
functions. 
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g 



Figure 1: Qualitative description of the Hubbard model spectrum. The upper line is the AF state. 
The lowest dashed lines in the perturbative multiplet are the light states. We also show two states 
from a different multiplet. One of them crosses a light state. At large g all lines are linear in g. At 
small g, A{g) = L + 0{g^) in the perturbative multiplet. 
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Figure 2: Spin correlation function for the lightest |FS) state of the Heisenberg chain (under the 
cyclic and S = constraints). 
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Figure 6: Spectrum flow for the Hubbard model at L = 8. 
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Figure 8: Asymptotic free fermion states for the perturbative multiplet at L = 8. As explained in 
the text, the fermions encircled by dashed ellipses indicate one of the various (spin) components of 
the actual state. 
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Figure 11: State |FS) at L = 12. Coupling dependence of Afs(5')- 
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Figure 12: State |FS) at L = 12. Fermion configuration at strong coupling. 
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